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ABSTRACT 

<n ; 

A standard version of a kinetic instability for the generation of Langmuir 
y—i '. waves by a beam of electrons is adapted to describe the analogous instability 

due to a beam of neutrinos. The interaction between a Langmuir wave and 
Q\ '. a neutrino is treated in the one-loop approximation to lowest order in an 

^ ! expansion in in the standard electroweak model. 

It is shown that this kinetic instability is far too weak to occur in a suggested 

OO . 

application to the reheating of the plasma behind a stalled shock in a type 
II supernova (SN). This theory is also used to test the validity of a previous 
^! analysis of a reactive neutrino beam instability and various shortcomings of 

this theory are noted. In particular, it is noted that relativistic plasma effects 
have a significant effect on the calculated growth rates, and that any theoretical 
description of neutrino-plasma interactions must be based directly on the 
electroweak theory. The basic scalings discussed in this paper suggest that a 
more complete investigation of neutrino-plasma processes should be undertaken 
to look for an efficient process capable of driving the stalled shock of a type II 
SN. 



Subject headings: elementary particles — instabilities — plasmas — stars: 
supernovae: general 



1. Introduction 

Processes involving neutrinos are important in a variety of astrophysical contexts, 
usually involving highly degenerate, relativistic plasmas. For example, neutrinos provide a 
cooling mechanism for red giants in the later stages of their evolution through the generation 
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of neutrino-antineutrino pairs from photons (Adams, Ruderman and Woo 1963; Braaten 
and Segel 1993). Of particular interest here is the intense flux of neutrinos generated in the 
first few seconds of a type II supernova (SN) explosion (Bethe 1990). There is a difficulty in 
understanding how the explosion occurs. The SN shock, which is needed to eject the stellar 
envelope, stalls due to energy losses through dissociation of nuclei. Heating by neutrinos is 
thought to be required to drive this shock. Wilson (1985) considered heating through direct 
neutrino scattering: highly beamed neutrinos are scattered by particles in the post shock 
plasma and deposit energy and momentum into this material, thereby driving the shock 
expansion. 

An alternative mechanism for the deposition of the neutrino energy behind the SN 
shock is the neutrino beam instability (Bingham et al. 1994), where the neutrinos emitted 
from the SN core drive a plasma instability which produces Langmuir turbulence. These 
Langmuir waves are then collisionally damped, thereby heating the plasma, and driving 
the SN shock. This instability may be understood by analogy with a more familiar 
electron beam instability: neutrinos propagating through a plasma emit longitudinal waves 
(Langmuir waves in an unmagnetized plasma) just as electrons do. This implies that for 
any plasma instability which relies on Cerenkov emission of Langmuir waves by electrons 
as the driving microscopic process, there is an analogous neutrino driven plasma instability. 
The simplest plasma instability is the beam instability, where a beam of electrons (e.g., 
Briggs, 1964), photons (Gedalin and Eichler 1993; Melrose 1994), or neutrinos (Bingham et 
al. 1994), causes Langmuir waves to grow through induced Cerenkov emission. 

In this paper a systematic development of the neutrino beam instability of Langmuir 
waves is developed based on the analogy with the electron beam and photon beam cases. 
A kinetic version of the neutrino beam instability is described, and the growth rate for this 
instability is calculated. The results are applied to the SN shock problem, and it is argued, 
contrary to the conclusions of Bingham et al. (1994). that the instability is ineffective and 
unimportant. 

It is known that a neutrino propagating through a medium acquires an induced charge 
(Oraevsky and Semikoz 1987; Nieves and Pal 1994) through resonant interactions with the 
electrons in the medium. It is through this induced charge that neutrinos couple to plasma 
waves. Emission by a neutrino has been considered both for longitudinal waves (Tsytovich 
1964; Oraevsky and Semikoz 1987) and for transverse waves (D'Olivo, Nieves and Pal 1990; 
Giunti, Kim and Lam 1991). (We note that the emission of a transverse wave by a neutrino 
is kinematically allowed only if the refractive index of the waves is greater than unity). The 
interaction between a neutrino and a Langmuir wave is treated here using the standard 
electroweak model in the one-loop approximation. The derivation of the probability for 
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Cerenkov emission by a neutrino of a Langmuir wave to lowest order in an expansion in 
1/Myy, where M w is the mass of the W boson, is summarized in Appendix A. We do not 
consider the effects of extensions to the standard model, such as allowing neutrinos to have 
inherent electric and magnetic moments or a nonzero intrinsic mass. 

Bingham et al. (1994; hereinafter BDSB) calculated the growth rate of a reactive 
version of the neutrino beam instability and predicted that the neutrino flux from the SN 
core would be sufficient to drive the instability. The theory used by BDSB was overly 
restrictive, in that it was assumed that the neutrino beam was monoenergetic and that it 
propagated in a cold plasma, two assumptions which increase the calculated growth rate. 
We also show that the "effective potential" neutrino-plasmon coupling used by BDSB does 
not incorporate all of the necessary electroweak physics, which may also have a dramatic 
effect on the growth rate of the reactive instability. We develop a kinetic version of the 
reactive instability proposed by BDSB, using methods explicitly based on the electroweak 
theory and drawing on techniques used in analysis of electron and photon beam instabilities. 

In section 2 we discuss the emission of longitudinal plasmons by neutrinos and relate 
this process to the theory of plasma instabilities. In section 3 the emission and absorption 
coefficients are evaluated for some specific neutrino distributions. It is shown that for a 
beamed distribution the absorption coefficient is negative over a range of values, implying 
that wave growth occurs for a sufficiently intense neutrino beam. In section 4, the 
absorption coefficient is evaluated for parameters relevant to type II SNe and it is shown 
that the kinetic instability is far too weak to operate. A comparison with the work of BDSB 
is then made, identifying important concerns for an explicit electroweak calculation of the 
reactive version of the neutrino beam instability. We also suggest that collisional damping 
completely swamps any tendency to wave growth. Natural units (c = h = 1) are used 
throughout, except where otherwise specified. 

2. Emission and absorption of Langmuir waves by neutrinos 

The treatment here of a neutrino beam instability is closely analogous with the 
earlier treatments of an electron beam instability (e.g., Melrose 1986, p. 53) and a photon 
beam instability (Melrose 1994). The time evolution of the occupation number of the 
Langmuir waves is described by a kinetic equation. In general, the kinetic equation is 
written in a semi-classical form involving the probability per unit time, iUM(k,p), that 
a particle (electron, photon or neutrino) emit a wave quantum in an arbitrary wave 
mode M. For present purposes, we restrict attention to the emission of the simplest 
longitudinal electromagnetic quantum in a plasma, a Langmuir wave. After a description 
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of the mechanism through which a neutrino in a plasma may emit a Langmuir wave, the 
probability of this process, which is derived in Appendix A, is written down for massless 
neutrinos and is simplified. 

2.1. The effective charge on the neutrino 

Cerenkov emission by a neutrino in a plasma occurs through resonant interactions 
of the type shown diagramatically in Figure 1. An electron in the plasma scatters off a 
neutrino via either W^-boson or Z-boson exchange and in doing so emits a wave quantum. 
This process is resonant when the incoming electron momentum is equal to its outgoing 
momentum. The probability of the emission of a wave quantum by a neutrino in a plasma 
is then calculated by averaging the probability of this resonant process over all the electrons 
in the plasma. 

Alternately, one may consider the processes shown in Figure 2, and interpret the 
internal electron lines as electron propagators averaged over the electrons and positrons in 
the plasma (this is the approach adopted in Appendix A). There are then two contributions 
from each of the diagrams of Figure 2. The first contribution is due to the vacuum 
contribution of virtual electron-positron pairs. This term may does not contribute here 
because the emission of a photon by a massless neutrino in vacuo is kinematically forbidden. 
The second contribution is due to the real electrons and positrons in the plasma, and is 
equivalent to the average over the resonant interactions shown in Figure 1. To lowest order 
in 1/Mw the two diagrams of figure 2 may be combined to form the 4-point interaction 
shown in figure 3. 

Nieves and Pal (1994) have calculated the induced charge on a neutrino in a medium. 
The matrix element for the interaction shown in Figure 3 (cf. Appendix A) may be written 
in the form of an electromagnetic interaction J /1 (/c)A /1 (/c), from which the 4-current, J M , for 
the neutrino is identified. The actual value of the effective charge on the neutrino, which 
follows from J°, is a fraction of order Gpn e /T e of the charge of the electron, where Gf is 
the Fermi constant, n e is the electron number density, and T e is the electron temperature 
as an energy. For plasma parameters relevant to the post shock region of a type II SN, this 
fraction is of order 10~ 13 . 

The only relevant kinematic restriction to the emission process shown in Figure 2 is 
conservation of 4-momentum. Denoting the 4-momenta of the neutrinos and the photons 
by (e, p) and (u;(k),k), respectively, conservation of momentum is satisfied (for |k| <C |p|) 
provided that the phase speed of the wave is less than the group velocity of the neutrino. 
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The phase speed of the Langmuir wave is W2,(k)/|k|, where ^(k) is its frequency, given 
by the solution of the dispersion relation of the plasma. A general form of the dispersion 
relation for Langmuir waves in a Fermi-Dirac plasma is written down in Appendix B below. 
For massless neutrinos the neutrino group velocity (to lowest order in l/M\y) is the speed 
of light. Induced effects lead to a energy- momentum relation for massless neutrinos given 
by (Nieves 1989) e — |p| + \f2Gpn e) where n e is the number density of electrons. Although 
this implies an effective mass and a phase speed less than that of light, the resulting 
corrections are negligible in the applications considered here. For example, for electron 
densities relevant to the post shock region of a type-II SN, n e = 10 36 m~ 3 , the induced mass 
is less than 10 _7 eV for a neutrino with energy ~ 1 eV. Thus, we may assume the simple 
vacuum energy-momentum relation for the neutrinos and the kinematic restriction on the 
emission of Langmuir waves is 



< V 9, 



where v g = 1 for massless neutrinos. 



2.2. Probability of emission 

We turn now to the calculation of the probability of emission of a Langmuir wave by 
massless neutrinos in a plasma. This expression is then simplified to lowest order in h in 
order to calculate the growth rates for the neutrino beam instability. 

The calculation of the probability per unit time of emission of a Langmuir wave by a 
massless neutrino is outlined in Appendix A (cf. Hardy and Melrose, 1996). The probability 
may be written, viz. equation ( |A12| ), 

w L (k, p) = ^ |k J £ 2 > (lkl)27r^ + co-e) [(2s - uf - |k| 2 ] (l - , (2) 

where Gp is the Fermi constant, a is the fine structure constant, Z(|k|) is a plasmon 
normalization factor (defined in Appendix B), with 

f 2sin 2 6W + \ for u e , 
CV \ 2sin 2 ^-i for^ - 



and where u = ujl(M) is understood. 



A convenient simplifying assumption in the calculation of the emission and absorption 
coefficients is to neglect the photon recoil. The recoil is the change in momentum, p, of the 
neutrino due to emission of a wave quantum, with momentum k, and the recoil is small for 
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|k| <C |p|, which also corresponds to w < e. Expanding the argument of the ^-function of 
equation (12) in |k|/|p| and retaining only the lowest order term gives 

«fc(k, P) = ^|k| 2 ^Z(|k|) (l - ^ 2 8{u - k • v), (4) 

where v = p/e. Ignoring the small induced mass of the neutrino implies |v| = 1. Equation 
(H) may also be written 

™l(Kp) = C(k)5(cosxo - cosx), (5) 

with 

2a 

with cosxo — ^/|k| and where \ is the angle between k and v. 

Note that, with the neglect of the photon recoil, the probability per unit time of 
emission is independent of the initial energy of the neutrino. 



2.3. Kinetic Equations 

In the semi-classical formalism used here, the distribution of plasmons is described 
by their occupation number, N L (k), and the neutrinos are assumed non-degenerate, with 
distribution function normalized according to 

n ^ = Swf f ^ (7) 

where n v is the number density of neutrinos. The time evolution of the plasmon distribution 
is governed by a kinetic equation, which includes those processes which change the number 
of plasmons present at a given energy through spontaneous and stimulated emission and by 
absorption. Thus, the kinetic equation may be written 



dN L (k) 
dt 



J d 3 p w L (k, p) [{1 + N(k)} /(p) - N(k)f(p - k)] . 



where the first term within the brackets represents spontaneous and stimulated emission, 
and the second term represents the inverse process, absorption of a plasmon (microscopic 
reversibility requires that the probability of emission and absorption be equal). 

For sufficiently smooth distribution functions, /(p — k) may be Taylor expanded, and 
then equation (||) gives 

^^ = a L (k)- 7L (k)iV L (k), (9) 
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where 



a L {k) = J d 3 pw L (p,k) /(p) 



(10) 



and 




are the emission and absorption coefficients, respectively. 



These coefficients appear in explicit solutions to equation (0). In the long time limit, 
the occupation number of the plasmons goes as 



N L (k) oc e 



7i(k)« 



(12) 



For stable distributions, such as a thermal isotropic distribution, 7z,(k) is positive, and 
the plasmon distribution tends to a constant level. If 7z(k) is negative, wave growth may 
occur. However, 7i(k) < is not a sufficient condition for wave growth, as there are 
other competing plasma processes which remove Langmuir waves from the plasma. Two 
of these are Landau damping by thermal electrons and collisional damping, with damping 
coefficients denoted by and 7 co n, respectively. For growth to occur, j L < must satisfy 



This condition sets the threshold on the neutrino intensity required for the instability to 
produce Langmuir turbulence. A particularly severe restriction is imposed by collisional 
damping which is of order (u p /A) In A, where A is the number of thermal electrons per 
Debye sphere. For the dense plasma around the SN core A is not a particularly large 
number (A ~ 10 2 ), so that collisional damping is a strong effect to be overcome in order for 
the Langmuir waves to grow. 

We turn now to evaluating the emission and absorption coefficients for a distribution 
of massless neutrinos. 



We assume that the neutrino distribution is axisymmetric and separable in spherical 
polar coordinates, (p, a, 0), where a = is the streaming direction: 



lL\ > 7coll + 7ld- 



(13) 



2.4. 



Emission coefficient 




(14) 



The plasmon momentum is written in terms of spherical coordinates (k, 8, <j)'). 
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To perform the integrals in equation (|T0| ) explicitly we need the following identity 
(Melrose & Stenhouse 1977), 

1 M \f**uxf \ f C ° sa+ Mcoso) 
d(cosa) / d(p b {cos xo- cos x) = -777 — ~ a r (15) 

-1 JO Jcosa- l<{a,t), XO) 

with 

F(a, 6, xo) = (1 + 2 cos a cos 8 cos xo — cos 2 6 — cos 2 a — cos 2 Xo) 1 ^ 2 , (16) 
and cosa± = cos(xo T Using equation (Pf), we find 

/•oo 

a L (k)=2C(k)g(e, X o) p 2 f(p)dp (17) 

JO 



with 



/•cos a+ 

g(8,Xo) = dcosa 

J cos a_ 



$(a) 



F(a,6>,xo) 



2.5. Absorption coefficient 

To evaluate the absorption coefficient, equation (|i~TD , we first write 



. d d k (cos a cos x — cos6>) d k 
k • — = k cos x~ — I z^TZ ~ l^"T 



dp dp p sin a da p sin a 

Our system is axisymmetric, so d/d(f> = 0. Thus equation fllTD gives 

7L (k) = -c(k) y o p 2 dpj^d 



d 

tttCosx 



£■ (19) 



cosa 



x 



J\(f)5(cosxo-cosx) (Ap^ + A9^jf(p,a). (20) 



The quantities Ap = kcosx and A9 = fc(cosacosx — cos#)/psina are the coefficients of 
the relevant derivatives in equation (|19|). Further evaluation of the absorption coefficient 
involves either performing the integrals directly or expanding the particle distribution in 
Legendre polynomials (Melrose & Stenhouse 1977). One finds 

7i (k) = 2|k|C(k) \2cos X o ~ sm X o^) g(d,Xo) f™ pf(p)dp. (21) 

Neutrinos of all energies interact with the waves satisfying x — Xo- F° r very low energy 
neutrinos, the condition |p| ^> |k| is not satisfied. The contribution from these low-energy 
neutrinos is unimportant in the present discussion, and the integral is extended to |p| = 0. 
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3. Specific neutrino distributions 

Equations ( |ITD and fl2T| ) may be evaluated for any separable axisymmetric particle 
distribution. Our main application of this work is to the neutrinos emitted from the core 
of type II SNe, where the neutrinos passing through a point in the post shock region are 
beamed due to their distance from the SN core. However, we first show that no isotropic 
distribution of massless neutrinos can cause kinetic wave growth. 



3.1. Isotropic distributions 

An isotropic distribution of neutrinos with arbitrary distribution function, f(p) has 
$(o) = 1. Substituting this into equation fllTf ) we find that 

g(e, X o) = rt, (22) 

and hence 

a L (k) = -C(k) / p 2 f(p)dp (23) 
I Jo 

and 

roc 

7i (k) = C(k)|k| cosxo / Pf(p)dp > 0. (24) 

Jo 

It follows that the absorption coefficient is non negative and hence wave growth is impossible 
for an isotropic distribution of massless neutrinos. This result also holds for electron and 
photon beam instabilities. 



3.2. Beamed distributions 

The analysis relevant to type II SNe is continued through the introduction of a 
particular form of beamed neutrino distribution which allows the analytic evaluation of 
equations (|17]) and (|2"1~1) . The absorption coefficient is then maximized (in the negative 
sense) with respect to the angle of emission to determine for which direction maximum 
growth may occur. 

We consider a highly beamed distribution of the form 

$(«) = -l exp j^rf, (25) 

where «o -C 1 is the beam opening angle (which is the angle subtended by the neutrinosphere 
on the sky for type II SNe). We need to evaluate g(9,Xo) f° r this distribution. Clearly, 
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g(9, Xo) is non-zero only for 9 — Xo ~ 0, so we expand F(9, Xo, ac) with a C 1 and 9 — xo ^ 1- 
This leads to 

F(9, Xo , a) « [a 2 -{9- Xo ) 2 ] 1/2 sin Xo . (26) 



Thus, 



and we also have 







d cos xo 



0(0. Xo) 



g(o,xo) 



271 



a sin Xo 



-V2^ 

3 ■ 2 1 

agsin xo 



-cxp 



<0 - Xo)' 



Xo)exp 



2al 



(27) 



(2? 



where we neglect the derivative of the slowly varying 1 / sin xo factor. Hence, the emission 



coefficient is given by equation (|T7|) with equation (|27D and the absorption coefficient is 
given by 

- Xo) 



lL {k) = 2C(k)|k| 2cosxo 



Q r 



9(9, Xo) / Pf(p)dp. 



(29) 



The right hand side of equation ( |29| ) is negative and largest in magnitude for 9 — xo ~ — a o- 
Thus the maximum growth for Langmuir waves occurs at an angle 9 w xo — a o to the beam 
axis. The maximum growth rate corresponds to 



7i m ax0) ~ _ 2g ( k ) fc s inxo^(xo - «o, Xo) / pf{p)dp. 



a 



(30) 



Equation (|30D and ( ^7|) imply 7i max oc I u /oiq, where I v is the neutrino intensity. 



4. Application to type II SN 

Neutrinos produced in the core of a type II SN are scattered many times as they 
diffuse outwards to a region which allows essentially free propagation (for a review of the 
physics of type II SN see Bethe 1990). Thus, the neutrinos have an approximately thermal 
energy spectrum and appear to be emitted from a surface enclosing the core known as the 
neutrinosphere. Hence, we may write the neutrino distribution function as 

f(p) = ^TP (31) 

With 120 

N = w T » ^ ( 32 ) 

which is defined by equation (|7|), and where T u and I v are the temperature (in units of 
energy) and intensity of the neutrinos, respectively. The use of a thermal distribution of 
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neutrinos implies that some neutrinos in the low energy tail do not satisfy the kinematic 
condition (|IJ), implying that the neglect of the photon recoil is not longer justified for 
these neutrinos. However, as already remarked, the fraction of neutrinos not satisfying 
this condition is negligible, provided the condition T v ^> Huj p is satisfied, as is the case in 
practice. 

At a distance r ^> r v from the core, where r v is the radius of the neutrinosphere, the 
neutrinos are confined to a ^ r u /r. Hence, we take 

r v 

a « — 

r 

in equation fl2"5"|). Evaluation of the integral in equation 



h^e~ l l 2 G 



2 „2 



7w(|k|)*^3 MVI w 




The Langmuir wave dispersion relation in the post shock region of a type II SN requires 
careful consideration. The relevant plasma is in an awkward parameter regime. Although 
degeneracy may be neglected, the plasma is neither ultrarelativistic nor nonrelativistic. 
Furthermore, the kinematic condition on wave emission, < 1, excludes the use of 
simple, known interpolation formulae which are valid only for > 1 (Braaten and Segel, 
1993). The approach adopted here is to calculate the dispersion relation, ui(ls.), and the 
normalization factor, Z(k), directly from the general expressions for an electron-positron 
plasma given by Hayes and Melrose (1984). This theory is summarised in Appendix B. 

We adopt the parameters used by BDSB for the material in the post shock region of a 
type II SN: n e = 10 36 m" 3 , I u = 3 x 10 33 Wm -2 , r = 300 km, and a = 3 x 10" 2 . We also 
assume a neutrino temperature of T u = 3 MeV. The electron temperature is thought to be 
around 2.5 x 10 9 K (Miller, Wilson, and Mayle, 1993). 



4.1. Absorption coefficient for SN plasma 

Using the techniques outlined in Appendix B, dispersion relations for a variety of 
plasma temperatures and a fixed electron density of n e = 10 36 m -3 are plotted in Figure 
4. This is a plot of phase speed, v^/c, against normalized wave number, |k|/2mc. The 
cutoff in the dispersion relation at high |k| occurs when the intrinsic damping of the mode 
(the Landau damping) becomes comparable to the frequency of the mode. In this event 
the plasma does not have a normal mode at that wave momentum, as any wave is highly 
damped within one wave cycle. In a nonrelativistic plasma, Landau damping limits the 
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range of phase velocities to > eV e = 3(T e /m) 1 / 2 . As T e increases, the range v m - m < < 1 
shrinks, where = t> m in is the limit imposed by Landau damping. In the ultrarelativistic 
limit the allowed range of < 1 becomes very narrow. 

This has immediate consequence for the strength of the emission process of Langmuir 
waves by neutrinos. Equation (FJ) shows that the emission probability is proportional to 
(1 — v 2 ) 2 which is small for ~ 1. This greatly suppresses the absorption coefficient in 
these high temperature plasmas. In the ultrarelativistic limit, the region for which < 1 
is vanishingly small, and all emission is forbidden. 

The maximum negative absorption coefficient, equation (|3~3|), for T = 2.5 x 10 9 K is 
plotted in figure 5 as a function of normalized wave number. It is clear that this growth 
rate is far too small to produce any Langmuir turbulence over the time scale of the 3 second 
neutrino burst. Indeed, this instability is too weak to overcome even the Landau damping 
of the SN plasma. The growth rate calculated here is approximately 11 orders of magnitude 
smaller than that calculated by BDSB, who obtained a growth rate for a reactive instability 
of ~ 10 6 s _1 . Part of the difference between these results lies in the the fact that a kinetic 
instability has been considered here (with a thermal distribution of neutrinos), rather than 
the reactive instability considered BDSB (who assume a monoenergetic distribution of 
neutrinos). The greatest part of the discrepancy lies in the supression due to the high phase 
speed of the longitudinal waves of the SN plasma, which is a factor of approximately 10 6 . 
This does not appear in the work of BDSB due to their artificial assumption concerning the 
properties of neutrinos, as is discussed below. 

One further point regarding the collisional damping of the plasma is required. This 
damping was not considered by BDSB. For the parameters of the post schock plasma 
assumed by BDSB the collisional damping is 7 co n ~ 10 15 s _1 , which is many orders of 
magnitude greater than the growth rates calculated either here or in BDSB. It is not 
clear to the present authors how any wave mode may be generated in a plasma with such 
overwhelming collisional damping. We conclude that the neutrino beam instability is 
ineffective in the post shock region of a type II SN. 

4.2. Comparison with BDSB 

The kinetic plasma instability discussed above is different from the reactive instability 
discussed by BDSB. As with most instabilities, these could be regarded as kinetic and 
reactive limits of a more general instability (e.g. Briggs 1964; Melrose 1986). The kinetic 
version applies when the growth rate is less than the bandwidth of the growing waves, when 
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the random-phase approximation applies, and the reactive version applies in the opposite 
limit when a fixed-phase disturbance grows in association with localized bunching of the 
particles (neutrinos). 

However, this equivalence does not apply to the kinetic instability discussed here and 
the reactive instability discussed by BDSB because of important differences in the underlying 
treatment of the neutrinos. The theory of BDSB describes the neutrino-Langmuir wave 
coupling through an effective potential in a Klein-Gordon equation for the neutrino wave 
field (Bethe 1986). Thus the assumed neutrino dynamics are those of a spin-0 particle, 
rather than the left handed spin- 1/2 particle of the standard model. The consequence of 
this assumption is that emission of a Langmuir wave in the forward direction is allowed 
in the theory of BDSB, whilst it is forbidden in the electroweak theory due to parity 
considerations. This is demonstrated in equation (01) above, and in D'Olivo, Nieves and Pal 
(1996). 

We commence our comparison of the theories from equation (20) of BDSB, which may 
be written in our notation as 



(ft + 27 ld -w = — - — p - / '—-d d p 35 

zna J w — k ■ v + zO 



where uj p is the plasma frequency, ft = uj — ry includes the real wave frequency and the wave 
damping, and iO denotes an infinitessimal imaginary part (the Landau prescription). This 
equation describes the contribution of neutrinos to the dispersive properties of Langmuir 



waves in the theory of BDSB. Note that in the limit of no neutrino flux, equation (pq) gives 
lu = uj p , which is the cold plasma dispersion relation, and 7 = 71a, which is the Landau 
damping of the finite temperature plasma. 

The growth rate of the reactive instability of BDSB is calculated by considering the 
real part of the right hand side of equation (p5|). The growth rate for the kinetic instability, 
in the cold plasma limit, is calculated from the imaginary part of the right hand side of 
equation (|35"D, due to the pole in the integrand. This imaginary part may be isoloated 
through the Plemelj formula 

V m5(uj — ijJq). (36) 



UJ — UJq + iO uj — UJq 

This imaginary part picks out the resonant interactions given by the condition uj — k • v = 

To illustrate the effect of the different assumptions made concerning the neutrino 
properties, we calculate the kinetic counterpart of the reactive instability considered by 
BDSB. The reactive version is obtained from equation ([35]) by omitting the factor iO 
and partially integrating before setting /(p) oc <5 3 (p), and solving the resulting algebraic 
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equation for solutions with complex frequency. The analogous kinetic version is obtained 
by retaining only the resonant part when equation (36) is inserted in equation (35). This 
gives 7 = 7id — 7^ DSB where 

7! DSB = / k • ^ ^ - k ■ v) d 3 p. (37) 



4a J dp 

The result (j37D is to be compared with equation (|Tl"|) in the cold plasma limit, where 
Z(k) = 1, and oj = u p , which gives 



GWv^M 2 L ^ V r, df(p) w3 rq 



For |k| pa cj p , as assumed by BDSB, equations (|37|) and (|38|) differ only in the factor (1 — v 2 ) 2 
in equation (^). This is precisely the term which appears through the absence of a right 
handed component ot the neutrino plane wave state, as discussed above. In particular, 
in a mildly relativistic plasma where the only relevant Langmuir waves have v$ < 1 and 
1 — v<j) <C 1, the factor (1 — vfy 2 in equation (38) implies that the growth rate for left-handed 
neutrinos is very much smaller than for the boson-like neutrinos assumed by BDSB. 



It is clear from the difference between equations fl37j) and ([38]) for |k| pa uo p that 
any analysis of neutrino-plasma interactions must be based on an explicit electroweak 
description of the microscopic processes occuring in the plasma. This has the added 
advantage of allowing the analysis to be performed at relativistic temperatures, high 
densities, and for broadband neutrino distributions, which are the major approximations 
in the work of BDSB. An electroweak theory description of the reactive neutrino beam 
instability is in progress. 



5. Conclusion 

In principle, a sufficiently intense, highly beamed distribution of neutrinos generates 
Langmuir turbulence in a manner analogous to an electron or photon beam. The kinetic 
version of this instability, which is studied here, is found to be too weak an instability 
to be relevant in any known terrestial or space plasma. The reactive counterpart of this 
instability was suggested by BDSB in connection with the reheating of the plasma behind 
the stalled shock of a type II SN. However, for simplicity, BDSB used a boson-like model for 
the neutrino rather than the electroweak theory which allows only left-handed neutrinos. 
From a comparison of the results of our work, based on the electroweak theory, with that of 
BDSB, it is clear that the two theories lead to quatitatively different results for Langmuir 
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waves with phase speeds close to the speed of light, cf. equations (|37| ) and fl38|). As the 
case v </, < 1 and 1 — v$ <S 1 is relevant to the SN problem, an electroweak theory of the 
reactive instability needs to be constructed, allowing for finite neutrino bandwidth, a finite 
temperature plasma, and correct inclusion of the neutrino spin effects. Such a theory is 
presently being developed. 

However, given the strong collisional damping of the SN plasma it is unlikely that even 
the reactive instability can be sufficiently strong to operate in the way suggested by BDSB. 
It may be that other plasma instabilities driven by an intense neutrino beam operate in 
such plasmas. A thorough investigation of the consequences of the large induced current 
generated by a strong flux of neutrinos in a dense plasma is required. Such an investigation 
should investigate the effect of an intrinsic neutrino mass, and other extensions to the 
standard model of particle physics, to determine whether there are any conditions under 
which neutrino plasma instabilities can be of astrophysical significance. 

SJH acknowledges valuable discussions with Dr. R. Bingham and Prof. J. Dawson, 
and thanks the Rutherford Appleton Laborotory for their hospitality during a recent visit. 



A. Calculation of the probability 

The theory of the neutrino-photon vertex, as shown in figure 3, has been derived to 
first order in (D'Olivo, Nieves and Pal 1989; Braaten and Segel 1993), where it is 

shown that the scattering matrix for longitudinal plasmon emission may be written 

M fi = %-^= *7 M (1 - %)var{k)e v (k), (Al) 

where a^ v {k) is the polarization tensor for the electron gas, e v (k) represents the polarization 
4 vector of the longitudinal Langmuir wave, a is the fine structure constant, and 

^ _ | 2 sin 2 9 W + | for u e , 

\ 2 sin 2 9w — \ for u T . 

As we consider coupling to longitudinal plasmons, only the longitudinal part of the 
linear response tensor contributes to equation ( |A1| ) (Tsytovich 1964; Braaten and Segel 
1993). Thus we may write (Melrose 1982; Nieves and Pal 1989) 



oT(k) ^a L (k)L» v (k,u), 



(A3) 
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where 



IT 



k 2 



k v ui i kv k 2 u»u v k"k u 
+ 



k 2 — (ku) 2 [ ku ku (ku) 2 k 2 
and, from the dispersion relation for the plasmons (Melrose 1986), 



a L {k) 



(kuf 
Ho 



10 

Ho 



(A4) 



(A5) 



The second form of the above expression is evaluated in the rest frame of the plasma, that 
is, the frame in which the 4-velocity of the plasma is it M = (1, 0) M . 

We choose the temporal gauge and write the polarization 4-vector for the plasmons as 



k v - (ku)u h 
[{ku) 2 



,21 1/2 



' Ikl 



(A6) 



Explicit evaluation of |M/j| through equations (|A3|) and ( |A"6| ) gives 



with 



and 



M, 



2 Ana 

ZpuPv - KPii - Puk^ + (pk)g^ u + ie a u/3ixk a p p 



(A7) 
(A8) 

v wv rwv ' (A9) 

with to = £Ji(k) and where p^ = (e(p), p) M and k^ = (^(k), k) M are the 4-momenta of the 
neutrino and the plasmon respectively. Conservation of 4-momentum implies 



i —Y(i ^ k 



p k = h*. 

F 2 



(A10) 



The probability per unit time of the emission of a Langmuir wave from a neutrino may 
be written 

w L (p, k) = \M fl \ 2 ^^ 2 ^'(P) + «L(k) - e(p)), (All) 

where e(p) and a»x(k) denote the neutrino and Langmuir wave dispersion relations, 
respectively. 



Substituting equations ( |A7| ) and ( |A1U| ) into equation ([All] ) one has 



Wi(p,k) 



G%4 |k| 2 «; 3 
167ra 



-2n8(e ! + to - e) (2e-u) 



LO 



(A12) 
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which is the probability per unit time of the emission of a Langmuir wave by a neutrino 
propagation in a medium. This expression is consistent with that derived by D'Olivo, 
Nieves and Pal (1996) who calculated the spontaneous emission rate of both longitudinal 
and transverse wave quanta. 



B. Longitudinal photon dispersion relations in a Fermi-Dirac plasma. 

The general dispersion relation for the longitudinal mode of an isotropic Fermi-Dirac 
plasma may be expressed to first order in the fine structure constant as an integral equation. 
This has been done by Hayes and Melrose (1984) and we reproduce the relevant details of 
their work here. 

The frequency of a longitudinal wave at a given wavenumber is given by the solution 
for uj = c<j L (k) of the dispersion equation 

uj 2 = -jj, a L (k) (Bl) 

where a L (k) is the longitudinal response function. This function is given by 



e 2 n P uj 2 



e muj r 1 



« L W = + T^T^ 2 ~ Ikl 2 )^) -mcoS^(k) + m 2 S^(k)}, (B2) 

m|k| z 27r^|k| d 4 

where 

r <i 3 p m . „ . . 

"« = 2 /(d?7"< £ > (B3 > 

is the proper electron number density, and 

n(e) = n~(e) +n + (e) (B4) 

is the sum of the electron, n~(e), and positron, n + (e) occupation numbers. The plasma 
dispersion functions of equation ( |B2| ) are given by 

de 



SW(k) = I —n(e)\nA 1 , (B5) 
J m 

S<U(k)= f ^n(e)£lnA 2 , (B6) 

S^ 2 \k) = ( ^n(£)£ 2 lnA!, (B7) 



with 



(e+-£ + m)(e + -e-ui){e+ + e- u)(e + + e + u) 

l\l — 7 77 T7 T7 r, (do) 

(e_ — e + oj)\£- — s — oj)(s- + e — oj)(e- + e + uj) 
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and 



(e+ — e + ou)(e + + e — — e — co>)(e_ + e + cu) 

(e_ — £ + c<j)(£:_ + e — l<j)(£: + — e — cu)(£ + + e + lo) ' 



e± = (e±2|p||k| 



(B9) 



(BIO) 



Various simplifications for the degenerate, nonrelativistic and ultrarelativistic limits 
exist for equation (|B2j ) reducing it to either a simple equation or an algebraic transcendental 
equation (Hayes and Melrose 1984; Braaten and Segel 1993). However, in its full generality, 
equation (|B^) must be solved numerically for specific electron-positron distributions. 



One further plasma quantity is required for the analysis of section 4, namely the wave 
normalization factor, Z(k), this is given by 



1 , A*o d L 

1 + {k) 



(Bll) 



and is twice the ratio of electric to total energy in the wave mode (Melrose 1986). 
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Fig. 1. — Feynman diagrams representing radiative decay of a neutrino mediated by a plasma 
electron. 




Fig. 3. — Feynman diagram for radiative decay of a neutrino to lowest order in 
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normalized wave numtJJ^2j77,c 

Fig. 4. — Plasmon dispersion relations for n e = 10 36 m~ 3 . Note that emission of Langmuir 
waves by neutrinos is only allowed for v<j, < c. The dispersion relation terminates at high 
wave number due to the destruction of the longitudinal mode by Landau damping. 
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Fig. 5. — Maximum negative absorption for n e 



10 



36 m 3 



3 x 10** Wm 



300 km, 



«o = 3 x 10~ 2 , T v = 3MeV and T e = 2.5 x 10 9 K. The range of wave numbers is that for 
which emission of a Langmuir wave is kinematically allowed, and for which a longitudinal 
mode exists in the plasma (see Figure 4). The absorption coefficient is far too small to 
produce Langmuir turbulence in the post shock plasma of a type II SN. 



